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Tli©  unstable  "burning  of  a  bipropellant  rocket;  combustion 
chamber  is  investigated  end  a  study  mads  of  the  requirements  for  an 
automatic  closed  loop  control  circuit  to  stabilise  the  motor. 

"  q  bipropellant  combustion  chamber  equations  developed  by 
Dr.  L.  Croceo'*'  are  utilised  as  tho  analytical  description  of  the 

rocket  motor  burning  phenomena.     Squ&tions  aimilar  to  those  developed 

(2) 

by  Dr.  H.  S.  Tsien   are  used  for  the  oxidissr  and  fuel  supply  systems 

and  the  te#o  closed  loop  stabilising  circuits. 

1!he  stability  or  Instability  of  the  systea  13  demonstrated  "a^ 
the  use  of  a  special  plotting  diagram  in  the  complex  plane  suggested 
by  if.  Satche  as  a  means  of  handling  systems  with  time  lag,  and 
developed  for  this  use  by  H.  S.  Tsien,  Oils  involves  separating  the 
transfer  function  into  two  parts*  In  the  complex  plane  the  first 
portion  of  the  transfer  function,  the  exponential  variable  containing 
the  time  lag,  plots  as  a  unit  circle  as  the  complex  variable  p  is 
made  to  take  a  contour  enclosing  the  positive  half  of  the  p-plane. 
If  tho  remaining  portion  of  the  transfer  function  intersects  thia 
unit  circle,  the  rocket  motor  can  be  unstable  for  large  reduced  time 
lag;  if  it  does  not  intersect  the  unit  circle,  the  system  is 
generally  stable,  although  the  roots  of  the  exponential  coefficient 
in  the  positive  half  of  the  complex  plane  faist  be  investigated,  -Shis 


ill 

latter  requj     t  can  bo  convenient \,     ecoti  !*  had  by  t    id  of 
yquist      .m. 

:ons  for  the  feedback  circuit      rir  ■        the 
oxidizer  and  fuel  transfer  function  requirements  are  determine  . 

o  cases  of  stable  combustion  and  t-zo  cases  of  unstable 
coabnstion  are  analyzed.  One  \m  case  ia  stabilised  by  the 

t.ion  of  a  feedback  circuit. 
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SYMBOLS  AITD  IHBTftTIOnS 

t  e      tine 

"T  «s   Inatantaneoua  value  of  the  time  lag 

n  =   pre83ure  exponent  of  preaaure  dependonoe  of  the 

proeeaaes  taking  place  during  the  tine  lag 

P  ss   instantaneous  pressure  in  the  conbustion  chamber 

also,  p        w  X   +  l^  rs   root  of  the  characteristic  equation 

with  the  induced  tine  aa  the  independent  variable 

p  =   preaaure  in  the  combustion  chamber  in  at  sady 

operation 

<f>  as       (-n  -  p)/p    =      fractional  variation  of  pressure 

in  the  combustion  chamber 

a      •      • 

ffi£,  s^t  me     «=   inatantaneoua  rate  of  injection,  burning,  and 

ejection  of  propellants 

£  ss   common  value  of  the  above  ratea  of  injection, 

burning  and  ejection  in  ateady  operation 

M-  aa   (ft^  •  a)/m  sc   fractional  variation  of  injection 

rate 

1%,  «s   inatantaneoua  ma33  of  gases  in  the  combustion 

chamber 

O-  =  l&gl  a  m      gaa  reaidence  tine  in  stea&y  operation 


Q.  3=  O  _  *  f'  =   total  reaidence  time  of  propollant9 

in  steady  operation 

$„  «=   absolute  temperature  of  combustion  gases 


VI 

A  =3  reduced  amplification  coefficient 

(^  a  reduced  ahgul  r  fre  juan 

r  =  Instantaneous  mixture  ratio 

r  «s  mixture  ratio  in  steady  ope  rati  o» 

K  a       coefficient  related  to  the  variation  of 

temperature  with  mixture  ratio 

K0  a      non-dinonsional  reservoir  capacity  for  the 

oxidizer  line 

Kg  a      non-din an 3 iunal  reservoir  capacity  for  the 

fuel  line 

H  a   (r  -  1)/  2  (r  +  1)  •-'-     coefficient  r  i   ating 

the  deviation  from  unity  of  tho  mixture  ratio 

i  a  subscript  indicating  injector  end 

e  m  subscript  indicating  ejector  end 

o  a  subscript  Indicating  oxidizer 

-#-  a  subscript  Indicating  fuel 

z  *=  t/Q      a   reduced  tine 

^  *=   ^/g   ■   reduced  time  lag 

C 

C^  a   constant  relating  the  mass  flow  rate  to  the 

pressure  in  the  feed  system  lines 

J  a   ^A/2ApA.  0   «   inertia  parameter  in  the  line 

a    p/3A  p  a   pressure  drop  parameter 

JLt  F^        a    ratio  of  functions  in  (d  /da)  describing  the 

response  of  the  feedback  circuits 

g  a    H  +  2K  •  Ke"*^  =   a  convenient  grouj     :  terms 


▼11 


A1,  B1,  D1,  L1,  'A1     =   commies  constants  used  In  the  trial 

solution 

O(p)  c      system  transfer  function 

Si(p)t  gp^  **       conponenta  of  G(:>)  used  in  the  Taien-^itche 

Diagram 

A,  C|  Df  Hi  3      «   convenient  system  constants  defined  In 

equations  (16)  through  (21) 

L,  M  «=   convenient  grouping  of  systen  coziat      id 

variables  defined  in  equations  (22)  and  ( 

$  =   2"0 (p)^  -  if  (, >)«  ■  fiui/p  *  &Q  +  s^  * 

a  convenient  symbol  for  the  feedback,  circuit 
transfer  function 


I.    r: :::::;  ujctioi: 

docket  motors  often  exhibit  rough  con'bustion  cipher  as 
chugging  (low  frequency  oscillation  boloitf  100  cycles  per  second) 
or  screoing  (high  frequancy  oscillations  several  hundreds  of  cycles 
per  seoond).     She  first  of  these  phenomena  has  been  investigate  1 
sufficiently  that  several  authors  have  developed  editions  describing 
analytically  the  oscillation  processes  in  rocket  ;r.otors*   '.     vfoose 
equations  which  are  utilized  in  this  paper  are  deduced  by  L.  Croeco 
under  the  folloying  assumptions: 

1.  !Ihe  processes  of  transforming  the  propellents  into  high 
temperature  gases  require  a  certain  tine  called  the  tine  lag,  <    ; 
It  Is  the  tine  after  ?*hich  the  liquids  nay  be  considered  suddenly 
transformed  into  hot  gases.     $he  rate  of  nixing  of  the  propellents 
depends  essentially  on  the  injection  system  used  and  reclrcul  ation 
effects,  but  probably  not  on  chamber  pressure.     G3he  conditions, 
droplet  size,   temperature,  vaporization  and  eventual  chemical  trans- 
formation of  the  propellents  after  raising  depend  considerably  on 
chamber  pressure. 

2.  The  time  lag  is  considered  the  same  for  all  particles  of  the 
propellents;  the  residence  tine  13  considered  the  sane  for  all 
particles,  I.e.,  all  particles  travel  from  the  injector  end  'inhere 
burning  occurs  to  the  exhaust  end  in  the  same  residence  time  interval. 


■  led 

to  bet 


; 


t 

f    (    )  dt1  ■  consent  (1) 


t  -f 
where  p  is  the  chanber  pressure  and  f(p)  is  the  tine  rate  of 
preparation  for  combustion. 

4.   Pressure  variations  are  sufficiently  slov  throughout  the 
co  .bustion  chamber  that  a  unifom  pressors  can  be  used,  p  -  i-(t)* 

5*        ature  of  the  gases  is  assumed  to  vary  with  tine  and  space 
but  not  ';lth  pressors  oscillations* 

■3*      l.ysis  of  the  stabilising  circuit  itself  and  its  oscillations 
af footing  the  rocket  notor  are  not  investigate  . 

L.  ~rocco    has  demonstrated  that  a  rocket  conbustion  chamber 
i9  Intrinsically  unstable,  i.e.,  the  burning  phenomena  itself  can 
generate  oscillations  even  without  variation  in  the  injection  rate* 
?his  occurs  because  any  oscillations  of  the  pressure  yill  cause  an 
oscillation  in  the  tins  lag  by  oqiv-tioa  (1).  IJhua  if  the  tine  lag  is 
oscillating  it  is  immediately  seen  that  -/hen  the  tine  la^  is 
decrees  ins.  the  horning  of  the  particles  that  /ere  injected  later  will 
catch  up  rith  the  homing  of  those  injected  earlier;  this  inoraa 
the  rate  of  burning  vith  rosoct  to  tho  average*  SIj     .y,  the 
o  >oslte  happens  for  an  increasing  time  lag.  Thus,  if  the  inci 
or  decreased  rate  of  burning  should  coincide  with  a  pressure 
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increase  or  decrease,  self-excitod  ovulations  can  exist  resulting 
in  unstable  burning. 

For  the  bipropellant  case  assumption  3  that  O  „    is  the 
residence  time  of  all  particles  means  that  each  particle  burns  at 
the  injector  end  and  travels  to  the  exhaust  end  carrying  with  it  the 
temperature  developed  at  the  combustion  instant*     2hus,  the 
temperature  at  the  eshsust,  SLq,  at  the  time  t  corresponds  to  tha 
mixture  ratio  of  the  propellents  injected  at  the  earlier  time 
t-9ta  t-  f  -  Q  _.     Ihis  can  be  put  in  terns  of  the  reduced 
time  by  dividing  by    0  -  giving,  i  -  <£  -  1.     Sierefore,  in  effect, 
we  have  two  reduced  time  lags,    S      and  1.0.     It  is  the  second  time 
lag  that  introduces  additional  complications  in  the  formulas  for 
the  bipropellant  case  by  appearing  as  an  exponential  term  in  g^Cp). 

She  problem  of  stabilizing  a  combustion  chamber  is  simplified 
by  the  use  of  the  Tsien-Satche  Diagram.     If  the  system  transfer 
function,  G(p),  is  separated  into  g^(p)  »  e"        and  g^p)  m  all 
remaining  terms  in  0(p),  then  as  the  complex  variable,  p  «  A  ■*  ia), 
it  made  to  take  the  contour  enclosing  the  positive  half  of  the 
p-plane,  stability  can  be  conveniently  determined  if  the  g2(p)  curve 
lies  completely  clear  of  the  unit  circle,  g^(p).    A  l^rquist  Diagram 
of  the  denominator  of  g*?(p)  also  is  required  to  determine  if  there 
are  any  poles  in  the  positive  half  plane. 

lor  an  unstable  case  in  which  g»(p)  does  intersect  the  unit 
circle,  g^(p),  stability  can  be  insured  by  designing  a  feedback 
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circult  transfer  function  yhich  will  cause  the  gp(p)  curve  to 

^ove  clear  of  the  g^(p)  curve,  and  still  meet  the  requirements 

of  the  ITyquist  Disgran  concerning  poles  in  the  positive  half  plana. 

For  a  bipropellant  rocket  the  temperature  of  the  combustion 
gases  is  a  function  of  the  local  nirture  ratio.  As  the  nixture 
ratio  operates  away  fron  the  steady-state  value,  r,  the  fractional 
variation  of  the  temperature  varies  as: 


2hia  defines  K  which  turns  out  to  be  an  interesting  parameter 
of  the  systen.  K  nay  be  either  positive  or  negative,  although, 
as  pointed  out  by  L.  Orocco,  K  is  generally  positive. 


ii.    wmxmms  of  sqjjATioHa 

Under  the  assumptions  given  above,  the  equation  of  the  combustion 
chamber  for  a  bipropellant  rocket  was  developed  by  L.  Crocco:  (See 
Figure  1  for  schematic  diagram  of  system). 


ii  +  H)/*^   +   U-HljtLffc  +  h(f-  4/s) 


(2) 


Following  H.  S.  7  si  en  and  considering  the  propellant  to  be  fed 
by  a  centrifugal  pump  run  at  constant  speed  the  following  equations 
for  the  bipropellant  feed  system  are  developed! 

Squations  (2),  (3)  and  (4)  specify  the  transient  behavior  of 
the  fuel  and  oxidizer  supply  systems  coupled  with  the  combustion 
chamber  when  the  reservoir  capacity  in  the  feed  line,  K0  and  Kp, 
are  specified.  Ihe  closed  loop  stabilization  circuit  can  be 
completed  by  requiring  Kq  and  Kf  to  be  a  function  of  <P     .  Thia 
means  that  the  reservoir  capacity  responds  to  pressures  from  the 


Mon  chamber     through,  a  suitable  amplifier  and  3orvo  system, 
?ms  Kq  and  Xf  are  assumed  to  satisfy  a  pair  of  integro-differenl 
equations, 

J-o  <  ?£  )  4»      =  Ko  (5) 

*<£>*       =I*  (6) 

where  F.  and  F0  correspond  to  a  ratio  of  functions  (generally 
polynomials)  in  d/ds. 

Trying  a  solution  of  the  formj 

<j>--A'*+*  *.-»'**'  M.f-D'^ 

(V) 

where  now  p  «   A     ■♦•  i  cO      is  complex,  aa  are  A»,  B1,  D1,  L1  and    * . 

eliminating  the  eoomon  factor  e^z»  equations  (2),    (3),    (4),    (5) 
and  (6)  represent  a  solution  of  the  system  if  the  five  homogeneous 
first  dja^ree  equations  below  are  simultaneously  satisfied: 

<1  -  n  *  p  ♦  ne^  )A«-(g  +  JJe*^  B»  ♦  (g  -  i)****  «  o  <R> 


P^  ♦    [l  +£<P0+&)  +  J0p]  B»   J-  [><P0+*)p  t-  J0P2  ]*>•  - 
PfA«  +  [l  f  i(Pf^)  ♦  Jjp  ]DH   ^(P^)p  +  JfP2  Jm«  «  0 


(9) 

(10) 
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*0(P)  A«  -  L»  e  0  (11) 

?f  (p)  A«  -  M«  *  0  (12) 

Since  they  are  horaogeneous  equations,  they  can  be  simultaneously 
satisfied  if,  and  only  if»  the  follot/inc  determinant  of  the 
coefficients  of  the  above  equations  vanishes: 


£<^;  o        o       -/ 


Ff(i>)  °         °       o       -  / 


where  g  ■=   H  ♦  2K  -  Ke~I>  has  been  used  for  convenience. 

Expanding  the  determinant  gives  equation  (14)  below;  if  p 
satisfies  this  equation,  our  assumed  solution,  equation  (7),  will 
be  a  possible  solution  of  the  system. 


-6- 


-P,(3-i?[KiCa+i)]j  +fe-*l»*H^l 


-k  [  + 


y\ 


z>t<+m+vi  -  JoPf(n)  +%z(rv] +  * 


le^«E^  + 


fe   ffoq>)  ($*-*)<  x.  T+j    -    6^    =  ° 


(/*-; 
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For  the  system  to  "bo  stable,  the  roots  p  of  equation  (14) 
should  not  have  real  parts.     'The  problem  of  stabilizing  the  rocket 
motor  thus  becomes  one  of  designing  the  closed  loop  transfer 
function,  TQ(p)  **■  Ff(p)»  so  that  thia  criterion  is  satisfied,     Tais 
condition  can  "be  investigated  most  easily  "by  use  of  the  Tsien-Satche 
diagram  and  a  riyqulst  lUagram. 

If  the  complex  variable  p  is  aade  to  take  a  contour  as  in 
Tigure  3  enclosing  the  positive  half  of  the  p-plane  from  -  «o  to  ■*•  cO 
and  around  a  circle  in  the  positive  direction  vith  infinite  radius, 
each  root  of  CJ(p)  «  0  (equation  14)  in  the  positive  half  vlane  will 
cause  G(p)  to  make  a  complete  revolution  about  the  origin.     Separating 
0(p)  into  Wo  parts,  £■  (p)  and  gp(p)t  by  dividing  by  the  coefficient 
of  the  exponential  tera,  e      •%  givesi 


<S<f»  -    %,(?)  -  Uf> 


k[l*X,JfP+  (V+S*f)p+  C-hR^J-hfz<p)Uff(pM\\  °5) 


where 


A    =       1  +£<Ff  *     )  (16) 
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5      1  +  i  <Pfl  +  i)  (17) 


<*    **0 


0  =  n  +<P0+Pf )    (J*  /*    +5    +  £u)  *   (P0-Pf )    (H  +  2K  «►  £  +    £  )  ♦ 

Vf  <*♦&  >  +    £&♦-&>  (18) 

»=  [(JffJ0)  (ttf£)  +  i*f0+J0*f)  (*♦?)•*•  <V<rVf>  <H*  ^J    <19) 

H=      (Pf -P0)   <K+£)  (20) 

Sr      -  K  (Jf  P0  -  J0  Pf)  (31) 

M=     f  (H  +  2K  -i)[(A-D   (Bp+  JoP2)*  Jf  <Bp2>  + 

(22) 

IH    £  (H  +  3K  -f  I)  [  (B  -  1)   (AP  +  dFfP3)  +  J0(^p2)  +  J^3]* 

-  JfcTP  [(B  -  1)   (Ap  +  Jfp2)  +  J0   (Ap2)  ♦  JfJoP3]? 

(23) 
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JTote  that  In  these  definitions,  A,  B,  C,  D,  E,  S,  aire  constants 
of  the  system,  but  that  L  and  II  include  the  variable  (;p). 

'?ith  this  arrangement  of  terms,  0(p)  becomes  a  vector  with 
vertex  on  g^(p)  and  the  storting  point  on  g2(p).  As  p  rmkes  the 
stated  contour  of  Figure  2,  g^(p)  becomes  a  unit  circle  and  gg(p) 
prescribes  various  contours  dependent  on  its  form. 

It  is  also  necessary  to  investigate  the  coefficient  of  e 
(this  is  the  denominator  of  go(p)  as  written  in  equation  (15)  to 
ascertain  if  it  has  any  poles  or  aeros  in  the  positive  half  plnne. 

If  the  coefficient  of  e"*  P  has  s  poles  sad  r  zeros  in  the 
right  half  p-plano,  then  as  p  traces  a  clockwise  contour  as  in 
H^ure  2,  the  coefficient  will  trace  r-s  clockwise  revolutions, 
Therefore,  for  stability  gp(p)  oust  make  r  counterclockwise  revolu- 
tions around  the  unit  circle.  Ifaus  even  though  the  primary  condition 
of  stability  Is  that  go(p)  lid  completely  clear  of  the  unit  circle, 
g,  (p),  it  must  also  satisfy  the  conditions  dictated  by  the  number 
of  poles  in  the  positive  half  plane  of  the  coefficient  of  the 
exponential  term,  e~    • 
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III.     STABILITY  OF  3IPB0PELLAH2  SXSfZKS 

WITHOUT  FTJBACK 


Ito  investigate  first  the  coupled  feed  systen  and  combustion 
chamber  without  the  feedback  circuit,   set  K0(p)  end  Xf  (?)  equal  to 
zero.     Cliis  reduces  equation  (15)  to: 

lP  '  i  f »  Jo  Jf  ^  +  (D  ♦  Sl-Nfi  +  O  ««*  J  J  (34) 

where  the  constants  A,  2,  C,  D,  S  and  E  are  defined  in  equations  (16) 
through  (21). 

^_^j;!  (stable  Case,  K«  0.2) 

Let  n  «  0.1;  *<  -  1.0;  r  «=  2.75;  E  -  0.233;  *g  =  4500°Kj 

dT  /dr  ■  650;  X  =  0.2;  J0  -  4-°?  Jf  *  °-85  po  =  0-5s 

Pf  «  1.0. 

Substituting  these  valuos  into  equations  (16)  through  (21)  gives: 

A  a  2.5;  B  ~     2.0$  G  *  1.65;  9  ■  1.08;  E  ■  0.2; 

3  ■  0.72. 

5!hu3t   equation  (24)  beconea: 


4    G(f)-^r~ 


(£^0.9X1. 5"-/-  o-2p)(Z-o  +  +.ojb) 


-1 


o— 


How  let  p  as  i&)  ,  where   *)    is  the  reduced  angular  frequency, 
Squat ion  (25)  reduces  to: 


(-3.Z.^V  &+*oS)l  -/f.4fu>2-  +  4..5- 

(26) 


f-(/o*  ±orfz<r"u)i^+l]LS-ho.7.*^Jl\ 


Substituting  values  for    <^>      gives  the  Tsien-Satche  Diagram  for 
gg(p)  plotted  in  Figure  3.     Checking  the  denominator  of  g  (;)  for 
poles  gives  the  llyquist  Diagram  jtlotted  in  figure  4.     The  numerator 
is  plotted  in  Figure  5. 

Proa  the  figures  it  is  evident  that  this  case  is  stable  since 
the  go(p)  curve  is  completely  clear  of  the  unit  circle  and  nakes  no 
loops  about  it.     Also,  the  denominator  in  the  Uyquist  ;">ia^rao  is 
proven  not  to  have  any  zeros  in  the  positive  half  p-plane,   since  it 
does  not  make  any  revolutions  about  the  origin, 

"■■•?•■".  \\\      (Stable  Case,  K  «  0.3). 

Let  n    m'   0.1  j      c\    m  1.0;     f    =     3.75;     H    ■     0.333;   *_    -     4500°" .; 

6 

d-7g/dr    ■    1000;     K    =     .3;     J0    m     4.0;     Jf    ■     0.8;     iQ    ■     0.5; 
Pf    =     1.0. 

Substituting  these  values  into  equations  (16)  through  (21)  gives: 
A    be     3.5;     B    »     2.0;     0    ■     1.50;     D    *     0.36;     B    ■     0.3;     3    ■  1.08. 
Letting  p    «     io3  and  substituting  these  v;lue3  into  equation  (24) 
gives: 
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Again,   substituting  values  for  ^    gives  the  Tsisa-Satcho 
Diagran  for  g«(p)  plotted  in  Figure  6;   the  ITyquist  Diagram  of  the 

a 

denojiiaator  is  plotted  in  Figure  7.     The  numerator  is  the  a  .  in 

3aae  I  ("Figure  5). 

Changing  K  in  this  second  case  has  resulted  in  a  radical  change 
in  the  Tsien-Satche  Diagram.     The  curve  &9(v)  now  nakes  two  counter- 
clockwise loops  about  the  unit  circle,  although  it  foes  not  intersect 
the  curve  &i(v).     Nevertheless,  this  is  still  a  stable  case  since 
the  denoninator  of  go(p)  nakes  two  clockwise  revolutions  about 
origin  in  the  Hycpiist  Diagram. 

QA3B  lilt     (Unstable  Case,     X    «     0.2) 

Let  n    «    0.5;    c<_    «     1.0;     r    »     2.75i     E    »     0.233;   Tg    -     4500°K; 

dTg/dr    «    650;     X    a    0.2;     J0    *=     4.0;     Jf    e     0.8;     P      ■     0.5; 

If     a      1.0. 

Substituting  these  values  into  equations  (16)  through  (31)  gives: 

-  2.5;  B  a  2.0;  0  m     2.65;  S  *  5.72;  H  ■  0.2; 
3  «  0.73. 

Letting  p  «  izO   and  substituting  those  values  into  equation  (24) 
gives: 
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6(*)=e*-J('3*"  +"■***  W32t°^2y _i    (.8) 


Substituting  values  for  «**    sand  plotting  gives  the  5?si©n-3atche 
Diagrsn  of  figure  8  for  g2(p)»     ^ke  behavior  of  the  denominator  of 
g^(p)  2.3  3ho^n  in  the  plot  of  Figure  9,     Tho  numerator  is  plotted  In 
Figure  10. 

Inspecting  the  figures  it  is  seen  that  the  g0(p)  curve  now 
intersects  the  unit  circle  and  thus  the  rocket  motor  can  "be  unstable 
for  large  reduced  tine  lags.     This  shift  to  instability  is  primarily 
the  result  of  the  change  in  the  value  of  n  from  0.1  to  0.5.     Thus, 
"n",  the  pressure  exponent  of  pressure  dependence  of  the  processes 
taking  place  during  the  time  lag,  T*    ,   is  shown  to  be  an  ipr>ort?nt 
sy3ten  parameter. 

GASB  IYt      (Unstable  Case,     K  »     0.3) 

Let  n    «     0.5;       <*-     =1.0;  r    *     2.75;     H    «     0.233;     Tg    e     45GC 

d?  /dr    =     1000;     K    «     0.3;  J0    *     4.0;     Jf    =     0.8;     PQ    *=     0.5; 

*f    «     1.0. 

Substituting  these  values  into  equations  (16)  through  (31)  givasi 

A    ■     2.5;     B    m     2.0;     C    -  3.50;     13    «     4.50;     H    «=     0.3; 

S    »     1.08. 
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Letting  p  ss  i  Cx>  and  substituting  these  values  into 
equation  (24)  gives: 

_u5cj    f(-lz^/o^)c  -  tiz^-HZ-tr ) 

Substituting  values  for     ^     and  plotting  the  values  of  g0(i>) 
gives  the  £sien-Satche  Diagram  of  Figure  11.     T'he  ITyquist  Diagram  of 
the  denominator  is  plotted  in  Figure  12  and  the  numerator,  which  is 
the  same  as  in  Case  III,  is  plotted  in  Figure  10, 

In  this  ease  changing  K  from  0.2  to  0.3  did  not  have  the 
spectacular  results  as  in  the  stable  case,  although  the  gq(p)  curve 
does  have  a  more  complicated  shape. 
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IT.     SIJABILIMIOB  01  A  BirRDFELLMTT  ST 
BT  A  F2KBBACK  CIROUIf 


Turning  to  the  problem  of  stabilising  an  unstable  burning 
rocket  motor,   it  is  obvious  that  a  feedback  transfer  function, 
P0(p)L  +  ?f(p)Xt  is  required  that  will  move  the  63  (p)  curve 
completely  clear  of  the  unit  circle  for  snail  values  of  ^    plus 
satisfying  the  requirements  of  the  Hyquist  Biagram  for  all  values  of 

^  •     Since  the  g2(p)  curve  can  only  interseot  the  unit  circle  in 

.  -i«-> 

the  region  of  small  values  of  w    ,  the  e  terns  in  the  denomi- 

nator of  g2(p)  —  (equation  15)     —  can  be  expanded  in  a  power  series 
of  (.<■'<**/     and  higher  order  terns  neglected.     Also,  define  a  feed 
systen  transfer  function,  F  (p)L  +  7f (p)M»  as  a  series  in  powers  of  p: 

y0(p)L  +  Ff<p)M^  5  (p)  «=  a«x/p  +  «o  +  *LP  *  *&?  <3°) 

Substituting  these  two  changes  into  0(p),   equation  (15),  and 
letting  p  ■  icJ      gives: 


^/^= £*k.  5     fee/  *  o-H)j/:^-f  j^fcoj>j  r  3+ j.  (<«>n       i 


How  clear  the  denominator  of  go(p)  of  imaginary  ter^s  by 
multiplying  both  numerator  and  denominator  by  the  conjugate  of  the 
denominator  giving: 
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*> 


®iis  i3  the  approximated  equation  for  the  §g(p)  torn  of  the 


system  transfer  function  and  holds  for  snail   <^    only.     !Pems  of 

3 

order  (1  <*»  )  and  higher  order  were  neglected  in  this  development. 

How  It  is  necessary  to  specify  the  shape  of  the  stabilized  curve, 
g2(p),  near  the  unit  circle  and  derive  the  values  of  a_^,  a^,  a^,  ag 
that  will  neet  these  requirements.  Let  this  stabilized  cuxve  be: 

«g(p)  m   V  f  X(i^)  +  I(i«>>)2  +  S(i<*>)3  (33) 

Equate  this  to  gg(p)  in  equation  (31)  and  clear  by  multiplying 
the  denominator  by  equation  (33),  transposing  the  numerator  terns 
to  the  left  hand  side  and  collecting  even  and  odd  powers  of  (icx)  ) 
together: 
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[w[*U{  +  |-s]-»-x^-^^s]-f-Y[c^^J+  Wdi  +  ta,+  YaQ  +  Z«i](<:^r+ 


[iA  -U^B  +  JjTo0-n)]M*--+  {yl[c+K]+  A6(/-»)  +  vVa0  + X  «-,  }  -  o      (34) 


~~~^  (35) 

Wa, -f  ^ae-^  Yd_f  ](coo)  +    y*A_,    -    o  (?£) 

Hote  that  in  this  dovolo- riont  terns  of  order  (ia>)     and  higher 
'.lera  neglected,     Thn  coefficient  of  each  (i")  torn  rsist  be  independ- 
ently aero  for  these  eqoatio&a  to  be  satisfied,     therefore,   setting 

the  coefficient  of  each  term,   i.e.,    (i  uJ  )"  t    (i  oJ  )  ,    (i  cJ  )  , 

p 
(i  oJ  )~"  to  aero  and  solving  for  the  unknowns,  a_^,  s^,  a^,  a^  gives! 

to«l  ■  ° 

^  *  3&l.x  *  -  AB(l-n)  -  W  [c  +  S^ 

Waj,  +  Xao  +  la.!  «  -  [jQA(l-n)  +  Jfi3  (1««i)+Ba]  -If  f  D-B  +  sj  -x[c  «■  ?J 

•Jag  +  &i  +  Tao  +  Za>i  ■  -  [ J0A  «*  J£&KTfJ0   (1-n)]   -wfnJ^J^  +  f*  -s]  * 

-  X  [d  -  B  +  s]    -  T  [  S  *  2} 

(35) 
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2hese  four  equations  can  be  solvad  sizault^eously  to  give 
explicit  relations  for  a.^,  o^,  a^t  a.?  In  terms  of  system  constants, 
.  thus  F0<p)l  *  F^(p)Wt  the  feedback  system  transfer  function,  is 
determined* 

Solving;  for  the  unknoi/n  a_^,  %t  a^,  a^  gives j 

•©*-  [  w  &-**!    -  (o+  a) 

•4  «  -[  J0A(l-n)  ♦  J^Bft-a)  +  Ba]w  -  (D-BhS)  +  (C  ■*•  H)   (  X  -^  )  + 
[/03(l-n)£] 

(3?) 

For  snail    ^    the  approximated  ggtp)  equation  (32)  can  be  used 
with  the  feedback  transfer  function  equation  (30)  to  plot  the 
Tsien— Satche  Diagram*     for  large  ^    the  general  equation  (15)  cai3t 
be  used  in  plotting  this  diagram;  also  the  Uymlst  HLagram  of  the 
denominator  of  go(p)  ^mst  be  investigated. 

To  illustrate  the  procedure  outlined  above,   consitler  the  unstable 
case  discussed  previously  where  K  «  0,2.     Let  the  required  stabilised 
curve  of  g2(p)  be  described  as  follows: 


S2<p)  «W+X(iw)i  Y(U>  )2  +  Z(lcJ)3  a  -  2  ■»•  8^"  +  8U)  .  (38) 

2his  arbitrary  selection  of  thD  gu(p)  curve  is  mdo  to  place 
the  stabilized  curve  in  the  approximate  position  of  the  known  stable 
curves  of  Cases  I  and  IX.     Ifous,  f  *  -  2;  X  a  -  8;  T  =  -  8;   2  -  G. 
Also,  A  a  2.5;  B  ■  2.0j  G  a  3.65;   D  «  5.7:3;   R  a  0.2;   S  a  0.72  (other 
values  as  Case  III). 

Substituting  these  values  into  equation  (37)  gives:  a_^  a  0; 
^s-  2.6j  a^  s=  -  5.84;  a^  a  -  1.06.     "Hth  these  values, 
5  a  -  2.6  +  1.06oi3  -  5.84  io3. 

Using  equation  (32)  and  the  above  values  the  fallowing  expression 
is  obtained  for  gu(p)i 


(      s  _  ,   +3.0  -/&*<**■+<(*.*<«»    ]  (39) 


Solving  this  equation  for  various  values  of  <*)   gives  the  belov 
table  of  results! 


c*A 

V/+  XM  -nT&lV  2&'«*** 

o 

—  Zoo 

to 

-  f.<f2  -  c  o.f3 

-/?L-C  O.Zo 

1 

~/si-  c  2zr 

-/£T  -±  Z-oo 

X 

3 

_/22  -  i  3.5S 

- /./z  -  Cz'^ 

f 

i 

-f-o.<?7  -c  £.7 

— . , 
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This  table  shows  that  as  *->  becomes  greater  than  1/3,  tli© 
values  of  g?(p)  b&VQ   too  great  an  error.  Hoover,  from  figure  13 
it  Is  apparent  that  the  curve  has  "been  stabilised  satisfactorily 
for  small  ^  • 

It  is  now  necessary  to  investigate  the  behavior  of  g0(p) 
using  the  exact  equation  (15)  which  includes  the  feedback  circuits 
substituting  the  values  determined  above  into  equation  (15)  results 
in: 

-2p   j-/3.2o>a-f-2.5"  -f  c  (-Jico'-f/oPcJ)         | 


Substituting  values  for  <^    and  plotting  g^p)  gives  the 
5?sien-3atche  Diagram  of  Figure  14.     She  Hyquist  Diagrsn  of  the 
denominator  is  plotted  in  Figure  15;  the  numerator  is  the  same  as 
Case  III  and  is  plotted  in  Figure  10. 

From  Figure  13  it  is  evident  that  the  stabilised  curve  is 
completely  clear  of  the  unit  circle  for  small  <+>   •     However,  for 
large  <°     the  addition  of  the  feedback  circuit  has  caused  the  go(p) 
curve  to  change  in  the  Tsien-Satche  Diagram  from  no  loops  about  the 
unit  circle  in  Gase  III  to  one  with  two  counterclockwise  revolutions 
about  gi(p)  in  the  stabilised  curve  of  Figure  14.     Nevertheless, 
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%he  systen  is  still  stable  since  the  "yquist  Biagran  shows  that 
the  denominator  of  &,(?)  makes  t*.ro  clockwise  revolutions  about  the 
origin, 

Eaus  an  unstable  burning  rocket  motor  has  been  stabilised  ty 
an  appropriate  closed  loorp  feedback  system.     The  circuit  chosen 
to  stabilize  Case  III  was  a  combination  of  three  terms  of  a 
differentiating  circuit  and  one  torta  of  an  integrating  circuit.     3ie 
coefficient  of  the  integrating  circuit  turned  out  to  be  sero 
indicating  this  type  circuit  need  not  be  used. 

In  order  to  show  the  nature  of  the  feedback  function,  X.  (  .), 
set  3^(p)  =  0  and  solve  for  F  (p)  in  equation  (30). 

Vp>  s  <a.l/P^ao+alP*a2P2)  /  L  -    A-l/P  +  Ao  4  V  +  V3  i) 

livaluating  I  by  substituting  the  values  of  constants  used 
previously,  approximating  e***0      by  a  power  series  in  (i  *•->  )  and 
neglecting  higher  order  terms  gives: 

X.  -  -  10.59  (1  *>  f  -  2.33  (i^  )  (43) 

Equation  (41)  can  now  be  evaluated  by  using  the  known  values  of  all 
terns,  resulting  ini 

Fo(p)  .  L-  1.06(l**>r  -  5.84(1^)  -  3. Sol  /     L  +  10.59(1^)^  2.33(i<^>)3 

2 
«  -  1.13/   (i^)  *  2.57  -  12(1^)  +  54.77(1^)  (43) 

wherefore,  to  our  order  of  approximation: 
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A.x  «  -  1.12 

JL,       *   +     ?.57 

(44) 
Ai    »  -  12 

Ae>    xs  +  54.77 

It  is  now  necessary  to  choose  a  relationship  for  J0(p)  that  will 
approach  sero  for  large  to    and  still  retain  the  desired  stabilising 
characteristics  for  snail   a)      ♦     !Ehus  the  following  expression*  which 
is  arbitrarily  chosen,  describes  a  circuit  that  will  stabilise  the 
eonbustion  chambers 

I0(p)  -   (-1.13/p)   (1  +  ElP  /  If  Bgp2  ♦  Bgp3)  (45) 

Ho  evaluate  3^,  &>,  a,  in  tanas  of  the  constants  in  equation  (44), 
equate  equation  (45)  to  (43)  giving: 

B^  »  -  2.39 

Bg  *  -  10.71  (46) 

B3  a  +  23.9 
Thus, 

*0(p)  ■  (-1.12/p)  (1  -  2.39p/  1  -  10.71p2  +  23.9p3)  (47) 

This  equation  for  the  closed  loop  circuit  gives  the  desired 
stabilisation  to  the  rocket  notor  for  low  frequency  pressure 
variations  in  the  burning  process  which  could  have  caused  combustion 
instability. 
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